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We theoretically investigate the quantum dynamics of the center of mass of trapped àtoms, 
whose internal degrees of freedom are driven in a A-shaped configuration with the làsers tuned at 
two-photon resonance. In the Lamb-Dicke regime, when the motional wave packet is well localized 
over the làser wavelenght, transient coherent population trapping occurs, cancelling transitions at 
the làser frequency. In this limit the motion can be efïiciently cooled to the ground state of the 
trapping potential. We derive an equation for the center-of-mass motion by adiabatically eliminating 
the internal degrees of freedom. This treatment provides the theoretical background of the scheme 
presented in [G. Morigi et al, Phys. Rev. Lett. 85, 4458 (2000)] and implemented in [C.F. Roos et al, 
^~>j ■ Phys. Rev. Lett. 85, 5547 (2000)]. We discuss the physical mechanisms determining the dynamics 

and identify new parameters regimes, where cooling is emcicnt. We discuss implementations of the 
£\J ■ scheme to cases where the trapping potential is not harmònic. 
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I. INTRODUCTION 



The progress in làser cooling of àtoms and ions has set the stage for coherent control of the dynamics of quan- 
tum mechanical systems [1]. By means of làser cooling, states of the center-of-mass motion of trapped àtoms with 
high purity have been prepared [1-4], allowing for instance for their coherent manipulation for quantum information 
çC) \ processing [5]. Nevertheless, there is a continuous interest for new and efíicient cooling methods, which solve ex- 
perimental difficulties and increase the efficiency of the process. In this context, a laser-cooling scheme for trapped 
àtoms has been recently proposed [6], that exploits the principies of Coherent Population Trapping (CPT) [7] and 
allows to achieve almost unit probability of occupation of the trapping-potential ground state [6,8]. This method 
has been demonstrated to be an alternative to sideband [2] and Raman-sideband cooling [3,4], routinely used for the 
preparation of very pure states of the center-of-mass motion of trapped àtoms and ions. Furthcr applications of this 
cooling method (now labeled as "EIT cooling") has been discussed in several publications [9,10]. 

The focus of this work is to discuss theoretically the physical principies on which this method is based, and 
particular ly the role of quantum coherence between atòmic states on the mechanical effects of light on trapped àtoms. 
Thus, in the first section we introduce the electrònic level scheme composed by two stable or metastable states coupled 
by làsers to a common excited state, the A configuration, and discuss in general CPT when the transitions are driven 
by counterpropagating làser beams (Doppler-sensitive case). Here, we observe that in presence of an external potential 
Qi confining the center-of-mass motion, (transient) CPT is obtained when the làsers are set at two-photon resonance 
and the wave packet is well localized over the làser wavelength (Lamb-Dicke regime). In the second section, starting 
. , from a general approach we develop the theoretical model, assuming that the atòmic center of mass is confined by an 
external potential in the Lamb-Dicke regime: That allows to adiabatically eliminate the internal degrees of freedom 
and derive an equation for the external degrees of freedom only [11]. We discuss this equation in detail when the 
potential is harmònic, and derive a set of rate equations for the occupation of the vibrational states. Thereby, we 
identify the parameters regime where cooling is effective. In some limits, these equations reduce to the ones used 
in [6,8,9]. Nevertheless, a result of this paper is the identification of the bàsic mechanism characterizing the dynamics, 
that allows us to determine new parameter regimes where cooling can be efHcient. We discuss the limit of validity of 
the equations derived, give alternative interpretations of the dynamics, and consider possible extensions of the method 
to cases, where the center of mass is confined by a potential, that is not necessarily harmònic and whose functional 
form may depend on the electrònic state. 

We remark that the laser-cooling dynamics of trapped àtoms, whose internal transitions are driven in a A configu- 
ration, have been investigated in several works, as for instance in [12-15]. These, however, focused on different cooling 
mechanisms. This work, together with [6], extends these previous analyses to other regimes, characterized by novel 
features of the center-of-mass dynamics, as we discuss below. 



II. THE DARK RESONANCE AND THE MOTION 



In this section, we first discuss the internal dynamics and steady state of an atom whose electrònic bound states 
are driven by làsers in a resulting A-configuration. We focus on the conditions for which CPT occurs. Then, we 
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consider the center-of-mass degrees of freedom and discuss under which conditions the features characterizing the 
bare intcrnal dynamics are preserved, when the motion is taken into account. The discussion in this section and 
throughout the paper is restricted to motion in one dimension, identified here with the í-axis. This allows a simplcr 
cxposition without loss of generality. 




FIG. 1. (a) Level Scheme: The sòlid arrows represent the làsers at Rabi freqüències Qi, Q2, that couple to the transitions 
|<7i) — > e), í/2) — ► |e), respectively, and are detuned of A from atòmic resonance. (b) Addition of a probe at Rabi frequency 
fip and detuning Ap, coupling \gi) — ► \e). 



A. The dark resonance 



An exemplary atòmic level configuration where the effects of quantum intcrfcrcncc manifest is the A-transition. 
It consists in two electrònic transitions, formed by two stable or metastable states that we label |<?i), | <7 2 } 7 which 
are coupled by làsers to the same excited state \e). For a closed transition, the atom stops to fluoresce when the 
states \g\) and \g 2 ) are resonantly coupled (two-photon resonance), as shown in Fig. l(a): The system evolves into 
the dark state, a stable atomic-states superposition which is decoupled from the excited state because of destructive 
interference between the excitation amplitudes. This phenomenon is called Coherent Population Trapping [7], and 
the àtoms are found in the coherence (dark state) 

|*D> = H( fi 2|ffl>-íïl|fl2», (1) 



where íí = ^/Qj-|-f2| and íli (íl 2 ) is the Rabi frequency of the làser coupling to the transition \gi) — ► |e) (\g 2 ) —> \ e ))- 
Here, without loss of generality we have assumed fií, íl 2 real. The dark state is accessed by spontaneous emission, 
unless the system has been initially prepared in it. Thus, the density matrix /?d = |^ , d)(Í'd| is the steady-state 
solution of the master equation for the atòmic density matrix p: dp/dt = £ p, where £ is the Liouvillian defined as 

C p=L·l·l,p]+Kp. (2) 

Here, H — H + Vo is the Hamilton operator, and its terms have the form (in the rotating wave approximation and 
in the frame rotating at the làser freqüències) 

H = -hA(\g 1 )(g 1 \ + \g 2 }(g 2 \) 1 (3) 

Vo = \ (Í2i |e> ( 9l | + ü 2 \e) (g 2 \ + H.c.) , (4) 

where A = ui\ — íJl,i = ^2 — w l.2 ar e the làser detunings, with the atòmic resonance freqüències uij of the transition 
\gj) — > |e) and the freqüències of the corresponding driving làser wlj (j — 1,2). The operator K is the Liouvillian 
describing spontaneous emission, 
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Kp=- 1 -[\e){e\p + p\e)(e\]+ ]T 7 J ·|ft><e|p|e>< ft |, (5) 

where 71, 72 are the rate of decay into \g\), \g2), respectively, and 71 + 72 = 7- It can be easily verified that the dark 
state is a dressed state of the system, i.e. an eigenstate of H. The other two dressed states read [16] 

=cos6»|e)+sin6i|^c), (6) 
|V>_) =sin0|e) -cos^c), (7) 

where 

VA 2 + g - A 

|V-c) = ^(Oi| 5l )+n 2 |52)), (9) 

and where we have introduced the state IV'c), orthogonal to |e) and \tpr>)- The states (6), (7) are at eigenfrequencies 
5w± = (A =F \/A 2 + SI 2 ) /2, and since they possess non-zero overlap with the excited state |e), they have a finitc 
decay rate and are populated in the transient dynamics. We denote their linewidths with 7+, 7_. The steady state is 
accessed at the slowest rate of decay and, for later convenience, we introduce T , the time scale corresponding to the 
inverse of this rate. 

The drcsscd-statc picture is a useful tool for interpreting the atòmic spectra in a pump-probc experiment, where, e.g., 
a weak probe at Rabi frequeney ííp (Op <C Í2i,í7 2 ) couples to the transition \gx) — * |e) as shown in Fig. l(b), whilc 
its frequeney is let sweep across the atòmic resonance. Figure 2 shows the spectrum of excitation as a function of the 
detuning of the probe Ap, for a certain choice of the làsers parameters. Here, one can observe that the component 
of the spectrum at Ap = A is zero, corresponding to the situation where the system is in the dark state I^d)- 
Moreover, the spectrum exhibits two resonances centered at Ap = 5u±, whose widths correspond approximately 
(when |A|,0 3> 7) to 7+, j_, respectively, and can be identified with the dressed states |V>-) [17]. Note that 

these resonance have not a Lorentzian shape: The spectrum shares in fact many similarities with a Fano profile [17]. 
Typical excitation spectra, measured with a singlc ion in a trap, are reported in [18,19]. 




-10 1 2 3 A 

FIG. 2. Excitation spectrum /(Ap) in arbitrary units as a function of the probe detuning Ap in unit of 7. Here, Q = 7, 
A = 2.57 ílp = 0.057. 



B. The motion 



We consider now the center-of-mass motion in presence of a conservative potential, of which for the moment the 
form is not specified. Given the mass of the atom m, the momentum p, the position x and the potential U(x), the 
mechanical Hamiltonian is 

H mcc = ^ + U{x). (10) 
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We denote with \ifi e ) the eigenvectors of H mcc at the eigenvalues e. The full dynamics are now described by the Master 
equation 

^p=^[H,p\ + ÍCp = Cp, (11) 

where p is the density matrix for the internal and external degrees of freedom and 

H = H + H mcc + V. (12) 
Here, V describes the coherent interaction of the atòmic dipole with the làsers, and has the form 

V =\ (n 1 é klXCOS ^\e)(g 1 \+n 2 é k2XCOS ^\e}(g 2 \+K.c.) , (13) 

where the làsers are travelling waves at wave vectors fci and k 2l propagating along the directions forming the angles 
<fii, <fi 2l respectively, with the x-axis. In (13) the spatial dependence is explicitly included, which couples to the 
external degrees of freedom of the ion, while the Rabi freqüències íli, 2 are assumed to be constant over the spatial 
region where the ion is localized. The liouvillian JC describes the incoherent scattering processes, whereby a photon 
is spontaneously emitted under an angle (fi with the axis of the motion. It has the form: 

ÍC~p= -l[\e){e\p + p\e){e\] (14) 



- Y, U [ dcosM^cosfc-Kel [e ifc '»«» \e)( 9j \, 



where N (cos (fi) is the probability distribution for the angles of photon emission respect to the motional axis. 

In this system, at a given instant of time perfect destructive interference between excitation amplitudes occurs for 

the state 

|*d) - ^ (íÏ2|si,*> -íïie^ 100 -* 1 -* 200 -*') 1 !^,*)) , (15) 

where exp(i(fci cos <fi\ — k 2 cos (fi 2 )x) is the displacement operator, acting on the external degrees of freedom, and is 
a state of the center-of-mass motion. The state (15) is stable -and thus a dark state- if it is eigenstate of Ho + H mcc . 
This is always true when the làsers are copropagating and k\ cos (fi\ = k 2 cos (fi 2 (or, for one-dimensional motion as in 
this case, when the direction of propagation of the làsers is orthogonal to the axis of the motion, cos (fix — cos <fi 2 =0): 
Then, the motional state factorizes out in (15). For k\COS(fii ^ k 2 cos(j) 2 , on the contrary, one must consider the 
particular form of the confining potcntial. For instance, for free àtoms (U(x) = const.) a perfect dark state exists for 
fci cos (fií — ~k 2 cos 02 = k and reads |^d) = (SI2 |<7i , —Kk) — Cli\g 2 , Tik)) /fi. This property has been used to prepare 
very cold atòmic samples [20] . In presence of a confining potential, on the other hand, there exists in general no state 
\^d) that is perfectly dark. Approximate dark states have been discussed in [15] for a 1D flat-bottom and for a 2D 
harmònic trap. 

Nevertheless, transient CPT can be observed in trapping potentials and in Doppler-sensitive confi guratio ns when 
the àtoms are in the Lamb-Dicke regime (LDR) , i.e. when the size of their motional wave packet \J (Ax 2 ) is much 
smallcr than the wavelength of the light, k\_ 2 ^ (Ax 2 ) <C 1. In this limit, a hierarchy of processes in the excitation of 
the center-of-mass wave packet is estabilished. At zero order in £ = k í 2 ^/ (Ax 2 ), the cffccts due to the spatial gradient 
of the light-atom potential are neglected: the àtoms behave like they were point-like, and the coherent transitions 
take place at the làser frequeney (carrier). Then, after the transient time To the àtoms have accessed the internal 
dark state [^d)- At first order in Ç, effccts due to the finite size of the motional wave-packet become manifest, and 
transitions between diffcrcnt motional states (sidebands transitions) oceur. On this longer time scale, that we denote 
with Tç, the atom is optically pumped out of the dark state into another state of the motion. In the Lamb-Dicke 
regime, the relation » T allows for a coarse-grained description of the dynamics, where the internal state of the 
atom is assumed to be always the dark state [^d)- 

These arguments suggest that for a trapped atom in the LDR some of the properties of the excitation spectrum 
discussed for the Doppler-free case may also be applicable to the Doppler-sensitive one. Here, the carrier transition 
is predominant, whereas transitions which change the state of the motion (sidebands transitions) are of higher order 
in the Lamb-Dicke parameter, and can be interpreted as transitions due to a probe (fíp) set at the corresponding 
frequeney in the bare atom, as illustrated for instance in Fig. l(b). In the next section, we show that this interpretation 
is theoretically justified. 
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III. THEORY 



Hcre, we derive the equations for the center-of-mass motion in the limit where the LDR applies and whcn the center- 
of-mass motion is confined by the same potential at all three electrònic levels. The procedure consists in adiabatically 
climinating the internal degrees of frecdom from the dynamical equation at second order in the parameter £, and 
it corresponds to analysing the coarse-grained evolution on the time interval At such that 3> At 3> T . The 
formalism we use has been first developed in [11] for a two-level transition driven by a running wave, and later applied 
to standing-wave drives and multilcvel transitions in [21]. In the following, we outline the fundamental steps that are 
most general to all treatments, and refer the reader to [11,21] for details (we have used the same notation as in [21] 
when possible). 



1. Lamb-Dicke limit 



In the Lamb-Dicke limit Ç <C 1, the operators exp(ifcjx) appearing in (13), (14) can be expanded in powers of Ç. At 
second order in this expansion Eq. (11) can be rewritten as 



d _ 
lP = 



£ +£i+£2 



dt 

where the Liouvillians Cj describe processes at the jth order in the Lamb-Dicke parameter, and are defined as: 



(16) 



Cop = C p + ^[H mec ,p\, (17) 
Cip=^[xV u p[, (18) 
C 2 p=^[x 2 V 2 ,p]+ÍC 2 p. (19) 



Here, V\, V 2 are the first and second order terms in the expansion of V and read 

K = | E k > cos< ^ (l e >^l - \9i)(e\) , (20) 

3=1,2 

^ = "7 E fc?cos 2 J ·íï J ·(|e)< ft -| + | ft ·)<e|). (21) 

.3=1,2 

The Liouvillian JC 2 has the form: 

ÍC 2 p = a ljtf\9j)(e\ fapx - x 2 p - px 2 ) \e)(gj\, (22) 

.3=1,2 

where a = j_ x d cos çW"(cos 4>) cos 2 <\>. 

At zero order in £, internal and external degrees of freedom are decoupled: The state pst, solution of Cop = 0, is 
not uniquely defined, and has the form pst — Pst <S> /"(O), where p$ t = po is the internal steady state and p(0) = 
TrjntjPoPÍO)} is the reduced density matrix, calculated from p at t — by tracing over the internal degrees of 
freedom (Tr; nt {}) and applying the projector Vo acting over the external degrees of freedom. The latter is defined as 
VoP = E e ~^2ij e ,ip' e IV'eKV'ÉKV'eIpIV'É)! wnere IV'e), W f ) are eigenstates of H mcc at e. In general, at zero order equation 
dtp = Cop admits an infinitc number of stable solutions. They can be expanded in the basis of (left) eigenvectors 
/5 e>e ' = pst <8> \ipe)(ipe' \ at the (imaginary) eigenvalues X e ^> = — i(e — e')/h of the Liouville operator Co, satisfying the 
secular equation Cop f _, f j = \ e ^'p e ^' (pst is eigenvector at A = 0). The eigenspaces at the eigenvalues A e , e < may be also 
infinitely degenerate, as it oceurs for instance in the harmònic oscillator. For £ thcse subspaces are coupled by C\, 
C 2 . At second-order perturbation theory in £, for QjÇlj <C min eie '^ e (|e — e'\) (i.e. when the spectum of Co is sufficiently 
spaced, to allow for non-degenerate perturbation theory), a closed equation for the dynamics in the subspace at À = 
can be derived. Denoting with To the projector onto this subspace, defined as VoP — Pst <8> T^intiTop}, this equation 
has the form [11]: 
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V0C2V0+ I dTV Cie CoT dVo 



p(t). (23) 



After substituting the explícit form of C\,Ci in the second term on the right-hand side of (23) and tracing over the 
internal degrees of freedom, we obtain: 

-/i = -V ^ J dr(Tr int {Vie £ °^ip St } [z, [Í(t), M ]] 

+ Tr int {V lC c ^{ Vl ,p st }} [x,fix(r)]). (24) 

Here, the matrix fi = Tri nt {'Po/5} is the reduced density matrix for the external degrees of freedom in the subspace at 
eigenvalue (at zero order) A = 0. The operator x(t) is here defined as x(r) = exp(— iH mec T /h)x exp(iH mcc T /h) . 

It is rcmarkable that the term V0C2V0 = 0. This result is explained by looking at the form of (19). Whcn tracing 
over the internal degrees of freedom, the first term of (19) gives rise to a contribution proportional to Tr; nt {V2/9st}: 
This term usually gives rise to a shift to the cigcnvalucs A e . e ', it represents a rcnormalization of the harmònic oscillator 
frequency due to the presence of the làser ficlds, and here it vanishes since there is no occupation of the excited state at 
steady state. The second term in (19) describes the diffusion arising from spontaneous emission into other mechanical 
states [21]. Again, since at steady state there is no excited-state occupation, it vanishes. Thus, the disapperance of 
VqCíPo is due to quantum interference at zero order in the Lamb-Dicke expansion. 

For a non-degenerate spectrum of eigenvalues e, the reduced matrix fi is diagonal, and the equation for a matrix 
element has the form 

^(VeMV'e} =ÇC £ ,,% £ , £Í ) HVeN^> + (VvMVv)] + H.C, (25) 

where the coefRcient S(u) e , e >) is the value of fluctuation spectrum of the operator V\ at the frequency uj e ^> = [e — e'] /%, 
and reads 

1 f°° 

S(«> e , e >) = 72 / d ^ Tr int {Vic CaT V lPst } ^' h . (26) 

n Jo 

The coefficient C €i€ i = |(^ e |x|-0 e ') | 2 weights the coupling between the center-of-mass states \ip e ) and \tp e t) due to the 
photon momentum at second order in the Lamb-Dicke expansion. The equations necessary for the derivation of the 
explicit form of (26) are reported in appendix A. Equation (26) shows that the rate for the transition \ij) t ) — > \ij) t i) is 
given by the value of the excitation spectrum for a probe, whose interaction with the atòmic transition is described by 
V\ and which is detuned from the pump by bj e ^i (sideband transition). Here, the form of the potential enters explicitly 
through the coefhcients C e , e <, and implicitly through the assumptions on the spectrum that have lead to (25). 



A. Harmònic oscillator 

We now let the potential be harmònic at frequency v, U(x) = imv 2 ï 2 , and introduce the annihilation and creation 
operators a and of a quantum of vibrational energy Tiv, such that x = xo^+a), p = ipo(a^ — a), with x a = y/h/2mv 
and po — yjTimv/2. The center-of-mass Hamiltonian reads 

H mcc = hv [a^a + . (27) 

Now, |V>e) = \n) and e = hv{n + 1/2), where n = 0, 1, ... is the numbcr of phonon excitations, and the mechanical 

energies are equidistantly spaced by %v. The coefficients C n , n ' — xl(nS n \ n -i + {n + l)# n ',n+i) 5 an d thus at first order 

in the Lamb-Dicke expansion the relevant transitions between motional states are the blue sideband \n) — > \n + 1) at 

frequency uj^ — v, and the red sideband \n) — > \n — 1) at frequency u>l + v. We define the Lamb-Dicke parameter 

r]j = kjXo, that fulfills the relation Çi^ = ??i,2-\/2(n) + 1, with (n) average number of phonon excitations. 

For the harmònic oscillator the equations derived in the previous section simplify notably: Equation (24) gets the 

form 

T-^/z = XqS(v) \—a)afi + a/xa^] 

+ xlS{-u) [-aaV + a)fia] + H.c, (28) 
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and for the probability /i„. n = (n|/i|n) of the system to be in the number state |n), Eq. (28) turns to a rate equation, 
whose form is well known in làser cooling of singlo ions [23] , 

^/Un,n = V 2 [( n + l)(-A-Mn+l,n+l - A+^ n ) 

+n(A 

— /J"n,n — l,n— 1 )]■ (29) 

In our case of a three-level atom, r] = t]i cos 0i — 772 cos 02, and 

1 (tt^V iv 2 



A± = 2Rc[S*(T^)l = t - J 7 r ± [ 5 • (30) 

Equation (29) has the same structure as the rate equation derived for sideband cooling in a two-level system. Hcre, 
however, the rates A± describe the sideband excitation including the effect of quantum interference between the 
atòmic transitions. Equation (29) allows for a steady state when A_ > A + , which is fulfilled when A < (blue 
detuning) and SI > 2v, or when A > (red detuning) and Cl < 2v. The value of the trap frequeney v = Cl/2 separates 
two regimes: for v < Cl/ 2 it is the narrow resonance that determines relevantly the center-of-mass dynamics, whereas 
for v > Cl/2 the sideband transitions are at the frequeney range of the broad resonance [22]. We remark that 
ï] = for k\ cos 0i = fc 2 cos 02, corresponding to the Doppler-free situation. Furthermore, the Lamb-Dicke parameter 
entering into the dynamics is the one determined by the làser wave vector. The Lamb-Dicke parameter connected to 
spontaneous emission events (i.e. recoils because of emission into other states of the motion) does not appear, since 
the diffusion term vanishes at second order in the Lamb-Dicke expansion. 

In the following we assume k\ cos0i 7^ k^ cos 02 and A < 0, Cl > 2v (A_ > A + ). Insight into the dynamics can be 
gained from the equation for the average number of phonon (n(t)) = X)^Lo n l· l n,n(t), that is derived from (29) and 
has the form [23] 

±(n) = -W(n)+ V 2 A + , (31) 
where W = r] 2 (A- — A + ) is the cooling rate. The steady state value (n)^ reads 

4[ny4-K«/-A)] 2 + 7 V 

Woo ~ 4HA|(f! 2 -4^) ' 1 ' 

and is minimum when Cl 2 = Ai>(i> — A). This relation corresponds to setting the a.c. Stark shift 6v + of the narrow 
resonance \^+) at the frequeney of the first red sideband, <W+ = A — v. For this value, (n)^ 1111 ' = (7/4|A|) 2 : 
Hence, low temperatures are achieved for làsers far detuned from atòmic resonance. This corresponds to an enhanced 
asymmetry of the excitation spectrum, as the one shown in Fig. 2, where the two resonances have very different 
widths. 

For Su) + = A — v the cooling rate seales as 

W m!íx - í7 2 (fiiíï 2 /íí) 2 /7. (33) 

Thus, fast cooling is achieved for large Rabi freqüències and when Cli — O2 - The ultimate limit to W is set by 
the parameters that ensure the validity of the perturbative treatment here applied: This is vàlid for rjjClj cos 9 7+ 
(j = 1,2), with 7+ ~ 7 cos 2 9 linewidth of the narrow resonance, corresponding in the bare atom to the situation where 
the probe (the sideband) does not saturate the transition to \^f+). At 5ui + = A — v, one has 7+ ~ ^ 7 /4VA 2 + Cl 2 , 
which sets the fastest rate at which emeient làser cooling can oceur, W m!íx ~ 7 +/2. 

It is remarkable that these results do not depend on the branching ratio 71/72- In fact, in this limit the branching 
ratio enters the problem only through T . Nevertheless, a too large branching ratio affects the time scale T at which 
the transient steady-state is reached. 

In Fig. 3 we test the validity of the adiabatic elimination procedure for various vàlues of the Lamb-Dicke parameter, 
by comparing the results predicted by (31) with a full numcrical simulation. The parameters are reported in the 
caption. Full agreement between the two results is found for rj — 0.02 (771 cos <p\ = —772 cos 02 = 0.01). It should be 
mentioned that in [6] full agreement has been found for 77 as large as 0.2. On the other hand, those results have been 
evaluated for the case Cl\ <C 2 , and the small value of Cli ensured the validity of the perturbative expansion. 
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FIG. 3. Plot of (a) {n)oo and (b) W as a function of rjo- Dashed line: rate equation result; 'o': numerical simulation, 500 
trajectòries with the Quantum Monte Cario method. Here, rji cos <f>\ = r]2 cos (j>2 = í?o (í? = 2r;o), v = 2 MHz, 7 = 20 MHz, 
fií = SI2 = 17 MHz, A = 70 MHz, 71/72 = 1. In (a) the result of the simulation in agreement with the rate equation prediction 
((n) = 0.005) is indicated by the arrow. In (b) the rate W is in units of 7/2. 



1. Discussion 



We have shown that, by properly choosing the làsers parameters, one can achieve almost unity ground-state 00 
cupation with this cooling method (EIT cooling). The state |^d)|0) is equivalent to the ground state in sideband 
cooling, sincc it is only off-resonantly (weakly) couplcd to other states, and it satisfies the criteria of an approximatc 
dark state as discussed in [15]. 

From (30) one recovers the rates of Eq. (4) in [6] in the limit fli <C Sl 2 - We have shown that the same dynamics are 
encountered in more general situations, that do not impose a specific relation between the two Rabi freqüències. From 
the technical point of view EIT cooling proves again to be more advantageous than Raman sideband cooling (see [6,9]). 
Such advantage is mainly twofold. On one hand, in EIT cooling both làsers cool the atom, and a decay into one or the 
other channel does not affcct the cfficicncy of the process, whilc in Raman sideband cooling a finitc branching ratio 
gives rise to heating [24]. Anothcr important feature of EIT cooling is the disappearance of the carrier absorption 
duc to quantum interference. This effect implies the suppression of diffusive processes: Since in the coarse-grained 
evolution the excited state is effectively empty, processes, where the atom is scattered into other motional states by 
spontaneous emission, disappear at second order in the Lamb-Dicke expansion. That implies an improved emeieney 
with respect to Raman-sideband cooling, where instead such processes are present, as already discussed in [6]. 

It is instructive to compare the dynamics in EIT cooling with the dynamics of a trappcd ion at the node of a 
standing wave, as studied for example in [21]. At the node of a standing wave the carrier absorption cancels, since 
here the value of the elèctric field is zero. Nevertheless, sideband absorption oceurs because of the finite size of the 
motional wave packet. In the case of a A- configuration driven by two travelling waves at two-photon resonance, 
the transient dark state (15) is a superposition of the states \gi) and \g 2 ) whose relative phase is a function of the 
coordinate x, so that the finite size of the wave packet allows sideband absorption also in this case. Nevertheless, in 
the LDR the gradient of the phase over the wavepacket is small, and the sideband transitions are excited on a longcr 
time seale. This can be illustrated when writing the atom-laser interaction (13) at the first order in the Lamb-Dicke 
expansion and in the form 

V « [|e)<* c | + i|e)<* D |fc»; + H.c] , (34) 

where we have made the simplifying assumptions fl-^ = Çl 2 = íl/y/2, fcicos^i = — fc 2 cos0 2 = k. Here, we see that 
the dark state is coupled to the excited state at first order in the Lamb-Dicke expansion, for effects arising from the 
finitc size of the motional wave packet. 

The atom dynamics during the coarse-grained evolution can be interpreted in terms of ficld gradients over the size 
of the wave-packet, that give rise to forces [25]. In this respect, one can say that this method uses the phase gradient 
of the dark state, due to the spatial gradient of the total field, for achieving cooling. In this context, we remark that 
the operator V\ in (26) is the gradient of the potential (13) at x = 0, i.e. at the center of trap. 
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Finally, we apply the results obtained for the harmònic oscillator to the case of a genèric potential U(x). Several 
conclusions drawn in this section are applicablc to the case describcd in Eq. (25), when the mechanical Hamiltonian 
has a discrete spectrum, and the minimum distance between two neighbouring energy level is sufficicntly large to allow 
for non-degenerate perturbation theory. Làser cooling is here achieved for the same parameters as for the harmònic 
oscillator. However, the narrow resonance enhances transitions in a finite range of freqüències (< P), and 6ui + must 
be properly tuned, e.g. to the average valuc of the red sideband transitions freqüències. The process will thus be 
cfhcient undcr the condition that, for each motional state, there is a sufhcient number of red sidebands inside this 
range, so that the rate of cooling for a given motional state is larger than the rate of heating. 

An interesting question is how the dynamics are affected when the external potential depcnds on the electrònic 
state, and thus when 

U(x) = LMaOIsiXsil + U e (x)\e)(e\ + U 2 (x)\g 2 )(g 2 \. 

We consider first the case Ui(x) — U 2 (x), whilc U e (x) is -say- constant, so that the center of mass of the excited 
atom is not spatially confined and the spectrum of iï m cc at the state |e) is a continuum. Assuming that for U\, U 2 
the Lamb-Dicke regime holds, then at two-photon resonance and during the transient dynamics the atom is optically 
pumped into the (transient) dark state (1). However, during To the center-of-mass wave packet changes, since each 
eigestate of U\(x) (U 2 (x)) may have non-zero overlap with several eigenstates of U e (x). This effect constitutes a 
diffusion mechanism, that lowers the cooling efficiency and, outside of some regimes, can make it even impossible. 
Formally, for U e ^ Ui,U 2 , the formalism applied in this section is not applicable, since one cannot separate the time 
scales characterizing the evolution of the infernal and external degrees of freedom. 

In the general case of three different confining potentials the presence of a dark state cannot be excluded: that however 
depends on the specific form of the functions Uj (x) . 

IV. CONCLUSIONS AND OUTLOOK 

We have presented a systematic investigation of the center-of-mass dynamics of a trapped ion, the infernal transitions 
of which are driven by làsers in a A-type configuration and set at two-photon resonance. Assuming that the center-of- 
mass wavepacket is well localized over the làser wavelength (Lamb-Dicke regime), we have adiabatically eliminated the 
infernal degrees of freedom from the equation of the center-of-mass dynamics, and obtained a set of rate equations for 
the occupation of the motional states. We have identified the parameter regimes where efficient ground-state cooling 
can be achieved. The derivation here presented provides the theoretical background for the equations in [6,8,9] and 
extends the parameter regime to cases which have not been previously considered. As also discussed in [6], we have 
shown that diffusive processes, encountered in cooling with two-level àtoms or with effective two-level systems (Raman 
sideband cooling), are suppressed because of quantum interference between the dipole transitions at zero order in the 
Lamb-Dicke expansion. Cooling takes place because of excitations due to the spatial gradient of the elèctric field over 
the width of the motional wave-packet, that are due to the finite size of the wave-packet itself and oceur at first order 
in the Lamb-Dicke expansion. The motion can be said to be coolcd by both làsers, while the branching ratio does not 
affect in general the efficiency of the process. 

Finally, we have discussed the possibility to observe these dynamics for other types of potentials, that may depend 
on the electrònic state. 

This work opens interesting prospeets in the manipulation of the quantum center-of-mass motion of àtoms by using 
quantum interference in driven multilevel transitions, that is subject of on-going investigations. 
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APPENDIX A: CÀLCUL ATION OF S(v) 

The term S(v) in (26) is the Laplace transform at w of the corrclation function G(t), defined as G(r) = 
Tri nt {Vi(r)Vi(0)pst}, where Vi(t) = V\C C ° T . This is evaluated applying the quantum regression theorem [16,26]. In 
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the following, we derive the equations that are essential for this calculation. For convenience, wc introduce the vector- 
operator à whose components are defined as: à\ — \gi){gi\, 02 = \gi)(g2\, 03 — \gi)(e\, (74 = |e)(<7i|, (75 = | «72) {e | , 
(i6 = | e) (92 1 , £7 = |<?2)(gi|, ò's = | <?i } {^2 1 - The mean value (àj) — Tr{àjp} obeys the equations d(àj)/dt = M(àj) + B, 
where M, B are a matrix and a column vector, respectively, and are defined through the equations 

8 n 

E M ^> = -7i((^i> + (*2>) - i-^ «*») - (*4>) , 

J=l 

8 n 

E^ifc) - -72((ai) + (à 2 )) - i-£ ((à 5 ) - (<r 6 >) , 



Y,M 7J {&j) = +i^(à 4 ) - i^(à 5 ), 

and = "fiój.i + 72^,2 + 1^(^,3 — ^,4) + i%-(^j,5 — àjfi), with j = 1, . . . , 8 and 5j.k the Kronecker-delta. According 
to this definition, the steady-state vector is now ost = M~ 1 B. 

Using this notation, we rewrite the operator V\ in (20) as V\ — ai(<74 — 0-3) +a 2 {à & — ér 5 ), with ctj = \hkj cos<t>jflj/2, 
j = 1,2. The Laplace transform S(v) is then the sum of the Laplace transforms Sj(u) of the individual tcrms 
9j( T ) = Tr{(7 :) (T)Vi(0)CTst}, such that S(v) = ai(s^(v) — s^{v)) + a2(ss(h') — S5(^)), whcrc Sj(v) are given by the 
equations 



8i {y) = ÇL jfe (TrWkVmpst} + Is fc Tr{^(0)p st }) , 
with L matrix, L = [\v — M] _1 . 
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